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Abstract 

We give a closed form for the correlation functions of ensembles of asymmetric real 
matrices in terms of the Pfaffian of an antisymmetric matrix formed from a 2 x 2 matrix 
kernel associated to the ensemble. We also derive closed forms for the matrix kernel 
and correlation functions for Ginibre's real ensemble. 

Difficulties arise in the study of ensembles of asymmetric random matrices which do not 
occur in the more commonly studied ensembles of random matrices (e.g. Hermitian, unitary, 
complex asymmetric, etc.). These difficulties stem from the basic fact that the eigenvalues 
of real matrices can either be real or members of complex conjugate pairs. Due to this 
additional complexity, our understanding of ensembles of real asymmetric matrices lags far 
behind our understanding of the Hermitian ensembles. In spite of these complications many 
important quantities such as a closed form for ensemble averages and correlation functions 
can be expressed using analogous structures for their Hermitian counterparts. 

Here we will show that the correlation functions for an ensemble of asymmetric real ma- 
trices can be expressed as the Pfaffian of a block matrix whose entries are expressed in terms 
of a 2 x 2 matrix kernel associated to the ensemble. We find much inspiration from Tracy 
and Widom's paper on correlation and cluster functions of Hermitian and related ensembles 
[19]. However, instead of using properties of the Fredholm determinant to calculate the 
correlation functions via the cluster functions, we use the notion of the Fredholm Pfaffian to 
determine the correlation functions directly. A Pfaffian analog of the Cauchy-Binet formula 
introduced by Rains [10] lies at the heart of our proof. For completeness we will include 
Rains' proof here. In place of the identities of de Bruijn [6] used by Tracy and Widom 
we will use an identity of the second author [16] to compute the correlation functions for 
ensembles of asymmetric matrices. Rains' Cauchy-Binet formula has applicability in a wider 
context than just the determination of the correlation function for ensembles of asymmetric 
matrices, and we will use it to give a simplified proof of the correlation function of Hermitian 
ensembles of random matrices when (3=1 and (3 = 4. 

Next we will apply our theory to Ginibre's ensemble of real asymmetric matrices. Ma- 
trices in this ensemble are square (say N x N), asymmetric with iid Gaussian entries. By 
scaling the eigenvalues by y/N, as N — > oo, the support of the eigenvalues converges to the 
unit disk. This scaling allows us to analyze the asymptotics of the matrix kernel in the bulk, 
and we will discover a new universality class of matrix models. This will yield a closed form 
for the limiting correlations of points in the bulk. Two different kernels emerge: one when 
expanding about points on the real axis and another when expanding about point in the 
(non-real) complex bulk. In the latter case we find that the limiting (large N) correlation 
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functions of Ginibre's real ensemble in the bulk off the real line are identical to the limiting 
correlation functions of Ginibre's complex ensemble in the bulk. 

Acknowledgments. We thank Percy Dcift for helpful early discussions and for introducing 
the authors to each other. The second author would also like to thank Brian Rider for 
many helpful discussions regarding Ginibre's real ensemble. Finally, we would like to thank 
Eric Rains for allowing us to include his proof of the the Pfaffian Cauchy-Binet formula 
(Appendix B). 

1 Point Processes on the Space of Eigenvalues of R iVx7V 

We start rather generally since the results in this manuscript can be used to describe not 
only the statistics of eigenvalues of ensembles of real matrices but also the statistics of 
roots of certain ensembles of real polynomials. We begin with random point processes on 
two-component systems. 

Let X be the set of finite multisets of the closed upper half plane H C C. An clement 
£ G X is called a configuration, and X is called the configuration space of H. Given a 
Borel set A of H, we define the function Na ■ X — ► Z-° by specifying that Na(£) be the 
cardinality (as a multiset) of (£r\A). We define £ to be the sigma-algebra on X generated 
by {JJa ■ A C H Borel}. A probability measure V defined on £ is called a random point 
process on X. 

For each pair of non-negative integers (L, M) we define Xl,m to be the subset of X 
consisting of those configurations £ which consist of exactly L real points and M points in 
the open upper half plane H. That is, 

Xl.m -^{CeX: N K (0 = L and = M}. 

Clearly, X LjM is measurable and X can be written as the disjoint union 

X = [J Xlm- 

L>0,M>0 

Given a point process T on X, we may define the measure T'l.ai on X by 7l,m{B) := 
3>(.B n Xl.m) f° r eacn B e £ . The measure ?l,m induces a measure on Xl^m (and we will 
also use the symbol Tl^m for this measure). 

Given a matrix Y e M. NxN there must be a pair of non-negative integers (L,M) with 
L + 2M = N such that, counting multiplicities, Y has L real eigenvalues and M non-real 
complex conjugate pairs of eigenvalues. By representing each pair of complex conjugate 
eigenvalues by its representative in H, we may identify all possible multisets of eigenvalues 
of matrices in ~R NxN with the disjoint union 

X N := U X L .M- (1.1) 

(L,M) 
L + 2M = N 

Similarly, we may identify all possible multisets of roots of degree real polynomials with 
this disjoint union. Thus, when studying the statistics of eigenvalues of ensembles of asym- 
metric real matrices (respectively of the roots of degree N real polynomials) we may restrict 
ourselves to random point processes on H which are supported on the disjoint union given 
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in (1.1). That is, the eigenvalue statistics of an ensemble of real matrices is determined 
by a set of finite measures "Pl,m on X L ^ M for each pair of non-negative integers (L,M) 
with L + 2M = N. In this situation we will say that 3 is a random point process on X N 
associated to the family of finite measures {Vl.m '■ L + 2M — N}. 

From here forward we will assume that L, M and N are non- negative integers such that 
L + 2M = N, and a sum indexed over (L, M) will be taken to be a sum over all pairs 
satisfying this condition. Moreover, when we refer to a point process we will always mean 
a point process on Xn- 

2 Point Processes on Associated to Weights 

In this section we will introduce an important class of point processes associated to Borel 
measures on C. We will be particularly interested in measures which are a sum of two 
mutually singular measures: one of which is absolutely continuous with respect to Lebesgue 
measure on R and one which is absolutely continuous with respect to Lebesgue measure on 
C. The corresponding densities with respect to Lebesgue measure will allow us to construct 
a weight function which uniquely determines the associated point process on Xn- Point 
processes of this sort arise in the study of asymmetric random matrices and the range of 
multiplicative functions on polynomials with real coefficients. 

It will be useful to distinguish non-real complex numbers and we set C* = C \ R. 

We start rather generally by constructing measures on Xl,m from measures on R L x C M . 
The benefit in doing this is that it allows us to express important quantities associated to 
point processes (averages, correlation functions, etc.) as rather pedestrian integrals over 
R L x C^. To each (a,/3) el l x we associate a configuration in X L , M given by 

{a,/3} := {ai, . . . ,a L ,0i, ■ ■ ■ ,3m}, where [3 m = {(3 m , f3 m } n H. 

A given configuration £ <E Xl,m may correspond to several vectors in M. L x and we will 
call {(a,/3) : {ol,(3\ = £} the set of configuration vectors of £. 

A function F on Xl.m induces a function on R L x C* f specified by (a.,/3) F{a,f3}, 
and given a measure vl,m on R L x there exists a unique measure 7l,m on Xl,m specified 
by demanding that 

/ F(Od'? L MO= n J, 9M / F{a,/3}dv L M<*,0), (2-1) 
Jx l ,m L\M\Z m J r l xC m 

for every S-measurable function F on X. The normalization constant L\M\2 M arises since 
a generic element £ € -X"l,m corresponds to L\M\2 M configuration vectors. By specifying 
measures vl,m on R L x C* f for all pairs (L, M) and normalizing so that the total measure 
of Xn is 1, we define a point process on Xn- 

A very important class of point processes arises when we demand that the various vl.m 
are all related to a single measure on C. Given a measure v\ on R and a measure on 
C*, we set v to be the measure v\ + v 2 on C. We will write vl for the product measure 
of v\ on M. L and v 2 m will be the product measure of v 2 on C M . By combining vl and 
i>2M with a certain Vandcrmonde determinant we will arrive at the desired measures on 
R L x C* f . Given a vector 7 e we define V{^) to be the N x N matrix whose n, n' entry 
is given by 7^ _1 . We will denote the determinant of ^(7) by A(7), and define the function 
A : R L x C M -> C by 

A(a,/3) := A(ai, . . .,a L ,(h,0i, ■ ■ -,0m, (3m)- 
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Using these definitions we set vlm to be the measure on K L x given by 

-(a,/3) = 2 M |A(a,/3)|, (2.2) 



d{v L x v2m) 

and we will write IP£ M for the measure on Xl^m given specified by vl,m as in (2.1). If 
IP£ M is finite for each pair (L, M), then we set 

yV ^7 E where Z" := £ 5>L m (^l,m). 

(L.M) (L.M) 

We will call IP" the point process on Xm associated to the weight measure v and Z v will be 
referred to as the 'partition junction of IP". 

We set Ai and A2 to be Lebesgue measure on R and C respectively and let A := Ai + A 2 . 
If there exists a function w : C — ► [0, 00) such that f = wA (by which we mean dv/dX = w), 
then we will call w the weight function of IP". In this situation we set IP™ M := !P£ M and 
define IP™ and Z™ analogously. If P" has weight function u» then 

, r L M 

av LM 



d(X L x A 2M ) 



K/3) = 2 M n^^) II M/3m)[|A(a,/3)|, (2.3) 

^ 1=1 m=l J 



and we define Qi^m : R L x C M to be the function given on the right hand side of (2.3). 
The collection {SIl,m ■ L + 2M = N} plays the role of the joint eigenvalue probability 
density function, and we will call SIl,m the L, M- partial joint eigenvalue probability density 
function of IP™. 



3 Examples of Point Processes Associated to Weights 

Point processes on Xjy associated to weights arise in a variety of contexts. It is often the 
case that the weight function w is invariant under complex conjugation. In this situation, 
there necessarily exists some function p : C — > [0, 00) such that 

r pin) if 7 GM; 
10(7) := I 

[ p{l)p{l) if7eC*. 

When it exists, we will cal p the root or eigenvalue function of IP™ (depending on whether 
IP™ models the roots of random polynomials or eigenvalues of random matrices). The 
root /eigenvalue function is often a more natural descriptor of IP™ than w. In this situa- 
tion we will write and for IP™ and Z w . 

3.1 Ginibre's Real Ensemble 

In 1965 J. Ginibre introduced three ensembles of random matrices whose entries were respec- 
tively chosen with Gaussian density from R, C and Hamilton's quaternions [10]. Ginibre's 
real ensemble is given by R NxN together with the probability measure r\ given by 

dn(Y) = (2 7 r)- N2 /"e-^ T ^d\ NxN (Y), 
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where Xnxn is Lebesgue measure on R" x . This ensemble has since been named Gi- 
nOE due to certain similarities with the Gaussian Orthogonal Ensemble. Among Ginibre's 
original goals was to produce a formula for the partial joint eigenvalue probability density 
functions of GinOE. He was only able to do this for the subset of matrices with all real 
eigenvalues. In the 1990s Lehmann and Sommers [13] and later Edelman [7] proved that 
the L, M-partial eigenvalue probability density function of GinOE is given by SIl,m for the 
eigenvalue function 

Gin( 7 ) := cxp(- 7 2 /2) yj erfc( V2|Im( 7 ) | ). 

Consequently, the investigation of the eigenvalue statistics of GinOE reduces to the study 
of ?[ Gin l. 



3.2 The Range of Mahler measure 

The Mahler measure of a polynomial 

N N 



f(x) = J2 a nX N - n =a l[(x- ln ) (3.1) 

n— n— 1 

is given by 

M/) =exp|^ log | /(e 2 - 9 ) = |oo| JJ max{l, | 7n |}. 

The second equality comes from Jensen's formula. Mahler measure arises in a number of 
contexts, i.e. ergodic theory, potential theory and Diophantine geometry and approximation 
(a good reference covering the many aspects of Mahler measure is [8]). One problem in 
the context of the geometry of numbers is to estimate the number of degree N integer 
polynomials with Mahler measure bounded by T > as T — > oo. Chern and Vaaler 
produced such an estimate in [5] using the general principal that the number of lattice 
points in a 'reasonable' domain in R JV+1 is roughly equal to the volume (Lebesgue measure) 
of the domain. That is, if we identify degree N polynomials with their vector of coefficients 
in R N+1 and use the approximation 

#{g(x) € Z[x] : deg<? = N,fx(g) < T} « X N+1 {g(x) e R[x] : degg - N, (i(g) < T} 

= T N+1 X N+1 {g(x) G R{x] : deg 5 = N,n(g) < 1}, 

then the main term in the asymptotic estimate Chern and Vaaler were interested in can be 
expressed in terms of the volume of the degree N star body of Mahler measure, 

Un - {g(x) € R[x] : degg = N,fi(g) < 1}. (3.2) 

The volume of this set is given by 

(3.3) 



N + 1 



where ^(7) = max{l, |7|} _Ar_1 [5]. Consequently, the volume of the star body which leads 
to an asymptotic estimate of interest in Diophantine geometry is essentially equals the 
partition function for the random point process . 
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3.3 The Range of Other Multiplicative Functions on Polynomials 

We can generalize Mahler measure by replacing the function 7 1— > max{l, \ with other 
functions of 7. Given a continuous function </> : C — > (0, 00) which satisfies the asymptotic 
formula, 

0(7) - | 7 | as I7I -> 00, (3.4) 
we define the function <!> : C[x] — > [0, 00) by 

AT N 

* : ao JJ(a: - In) i-» l«0 1 J| 0(7n)- 
n— 1 n— 1 

The function $ is known as a multiplicative distance function (so named because it is a 
distance functions in the sense of the geometry of numbers on finite dimensional subspaces 
of C[x]). The asymptotic condition in (3.4) ensures that $ is continuous on finite dimensional 
subspaces of C[x] (one of the axioms of a distance functions). 
We define the degree N starbody of <j> in analogy with (3.2): 

Ujv($) = {g{x) e R[x] : deg 5 = N, $(<?) < 1}. 

In this situation the volume of LCjv( < &) equals 

\ N+1 (U N {*)) = j£ r[ zW, where = M) - " -1 - 

We may discover Further information about the range of values $ takes on degree N 
real polynomials by considering the point process on Xn corresponding to the root function 
1/1(7) = (j>{ , y)~ a for a > N. The partition function of CP^ is therefore a function of a, 
and ZW(ff) is known as the degree N moment function of $. In fact, we may extend the 
domain of a real moment functions to a function of a complex variable s on the half-plane 
Re(s) > N. Moreover, in this domain Z^(s) is analytic. Any analytic continuation beyond 
this half-plane gives information about the range of values of <I> which may not be realizable 
by other methods. For instance, when 1/1(7) = max{l, |7|}~ s , Chern and Vaaler discovered 
that Z^' (s) has an analytic continuation to a rational function of s with rational coefficients 
and poles at positive integers < N [5] . Similar results have been found for moment functions 
of other multiplicative distance functions; see [15]. 



4 Correlation Measures and Functions 

Suppose t and m are non-negative integers, not both equal to 0. Then, given a function 
/:K'xC™^Cwe define the function F f : X N -» C by 

{x,z}C£ 

where the sum is over all (x, z) € M. 1 x C™ such that {x, z} C £. We take an empty sum to 
equal 0, and thus if £ e X L>M with L < £ or M < m, then F/(£) = 0. 

Given a point process T on Xat, if there exists a measure on R e x C™ such that for 
all Borel measurable functions /, 



/ /(x,z)^, m (x,z)= / FfiOainO, 

JR e xC™ JX N 



(4.1) 
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then we call pe, m the (£, m)— correlation measure of IP. Furthermore, if pi tTn has a density 
with respect to X e x X 2m then we will call this density the £, m- correlation function of 7 
and denote it by Re, m - By convention we will take i? ,o to be the contant 1. 

Proposition 1. dpe, m — ( x z ) e q Ua l s 
d(vt x v 2m ) 



2 m x 1 



- J j |A(xVa,zV/3)|di/L_^(a)di/2(M-m)(/3), 



Z" ^ (L-£)!(M-m) 

w/iere xVael 1 is the vector formed by concatenating the vectors x e R e and a e R L ~ e 
(and similarly for zV/3 G C M J. 

Corollary 2. If v = wX, then ify m (x, z) equals 

E (£-l)!(M-m)!2^— / / ^,M(xVa,zV/3)dA L _,(a)dA 2(M _ m) (/3). 

Proof of Proposition 1. Assume that L > I and M > m, from (2.1) and (2.2), 

/ m)dnMO = TTm I I FM,P}\^{<*,P)\dM<x)dv2M{P), (4.2) 

The function (a,/3) i— > |A(a,/3)| is invariant under any permutation of the coordinates of a 
and /3, since such a permutation merely permutes the columns of the Vandermonde matrix. 
Similarly, replacing any of the f3 with their complex conjugates merely transposes pairs of 
columns of the Vandermonde matrix. That is, if (a, (3) and (a', 0) are elements in R L x 
such that {a, (3} = {a',f3'}, then |A(a,/3)| = |A(a',/3')|. Moreover, if any of the /3 are 
real, then the Vandermonde matrix has two identical columns and is therefore zero. We 
may thus replace the domain of integration on the right hand side of (4.2) with R L x C M . 
In fact, we may assume that the domain of integration on the right hand side is over the 
subset of R L x C M consisting of those vectors with distinct coordinates. 
Now, 

F f {a,/3}= £ /(x,z). 

{x,z}C{a„8} 

Assuming that the coordinates of (a,/3) are distinct, and if (x, z) G R e x C m is such that 
{x, z} C {a, (3} we may find a vector (a, b) el"x £M-m suc ^ ^at (x V a, z V b) is given 
by permuting the coordinates of (a,f3). Clearly |A(a,/3)| = |A(xVa,zVb)| and 

dv L {a.) dv2M{(3) = dvi{-x) dv L -i(a) dv rn (z) dv M -m(b). 

These observations together with an application of Fubini's Theorem imply that (4.2) can 
be written as 



\f\ J t j X] /( x ' z )| J j |A(xVa,zVb)| dv L - t {a)dvM- m (b) > du e (x) dv m (z). 



< K L-e C M 

Now, it is easily seen that there are 

' M 



in 
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vectors corresponding to each {x, z}, and thus we find 

JXl.m «* JC m 



1 {L-i)\{M-m)\ 
The proposition follows since 



|A(xVa,zVb)| dv L -t{a) dv M - m (b) > dv e (x) dv m (z). 



JX (L,M) L - M 

From here forward, and unless otherwise stated, (£, m) will represent an ordered pair of 
non-negative integers such that I + 2m < N. 



5 A Matrix Kernel for Point Processes on Xn 

From here forward we will assume that TV is even. 

Let "? v be the point process on Xjy associated to the weight w, and as before let v 
v\ + v 2 . We define the operator e v on the Hilbert space L 2 {v) by 



- g(y)sgn(y-j)du 1 (y) if 7 e E, 
z Jr 

ig(j) sgn(Im(7)) if 7 G C*, 

and we use this to define the skew-symmetric bilinear form on L 2 (v) given by 

{g\h)u-= / 0(7)^(7) - ^5(7)^(7)^(7)- 



JC 

If v = wX then (g\h) v — (g\h)\ where, for instance, ^(7) := h^w^). 

Theorem 3. Let q = {90(7), 92(7), • • • , QN-iil)} C R[y] be such that each q n is monic and 
degq n = n. Then, 

Z" = PfU£, (5.1) 

the Pfaffian of UjJ, where 

:= [{q n \qM; n,n' = 0,1, ... ,N — 1. 

We will call q a complete family of monic polynomials. 

Remark. It is at this point that it is necessary that N be even, since the Pfaffian is only 
defined for antisymmetric matrices with an even number of rows and columns. A similar 
formula to (5.1) exists for Z v in the case when iV is odd; see [16]. However, we have not 
pursued the subsequent analysis necessary to recover the correlation functions in this case. 
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Theorem 3 follows from results proved in [16]. In fact, [16] gives a formula for the average 
of a multiplicative class function over the point process on X N determined by Ginibre's 
weight. However, the combinatorics necessary to arrive at such averages is independent of 
any specific feature of Gin(7) and the measure on C specified by Gin(7)cZA(7) can formally 
be replaced by any measure v. 

In order to express the correlation functions for the point process V w associated to the 
weight w we will define 77 to be a measure on C given by a linear combination of point 
masses, and then use the definition of the partition function and properties of Pfaffians to 
expand both sides of the equation 



r/w(X+-q) 1 

Z = _Lpfutu(A+u) 



(5.2) 



The coefficients in the linear combination defining r\ appear again in terms on both sides of 
the expanded equation, and after identifying like coefficients on both sides of the expanded 
equation we will be able to read off a closed form for the correlation functions in terms of 
the Pfaffian of a matrix whose entries depend on a 2 x 2 matrix kernel. 

In order to define the 2x2 matrix kernel for CP™ , we let q be a complete family of monic 
polynomials, and we define 



Qn{l) ■= qn(l)w(j) n = 0, 1, . . . ,N - 1. 



(5.3) 



We then define 



N-l N-l 

Sjv(7,7') :=2^ ^n, 

n=0 n'=0 



yn (7) z\qn'{i), 



where we define n n ,n' to be the n, n' entry of (U™ ) T . Similarly we define, 



N-l N-l 



ISnM) :=2^ ^ 

(7) eA9n'(7') 



n=0 n'=0 



and 



N-l N-l 

ds n (~i^') -.= 2 Y X! Vn, n > q n {i)q n <(i)- 

n=0 n'=0 

Remark. The functions Sn,ISn and DSn can be shown to be independent of the family 
q. By setting q to be skew-orthogonal with respect to the bilinear form (-|-)a we arrive at 
particularly simple representations for these expressions. 

Finally, in order to define the matrix kernel f w we define the function £ : C 2 

by 

• 1 

if 7, 7' e 



{-|,o,±} 



£(7,7') := 



- sgn( 7 - 7') 







otherwise. 



The matrix kernel of J"" is then given by 



-s N {i, 7 ) is N { 1 ,i) + £( 7 , V) 



(5.4) 



Remark. The explicit A^-dependence of Kn and its constituents is traditional, since one is 
often interested in the N — > 00 asymptotics of K^. 
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6 Correlation Functions in Terms of the Matrix Kernel 

Finally we may state the main result of this manuscript. 
Theorem 4. The i, m-correlation function of O 31 " is given by 



ifc, m (x,z) = Pf 



K N (x j ,x jl ) K N (xj,z k >) 
K N (zk,Xji) K N (z kl z k >) 



3,3' = 1,2,..., £; 
k, k' = 1, 2, . . . , m. 



Remark. The matrix on the right hand side of this expression is composed of 2 x 2 blocks, 
so that, for instance, the first row of 2 x 2 blocks is given by 



[K N (xi,Xi) 



K N (xi,xt) K N (xi,Zi) 



K N (xi, z m )] . 



We define the measure S on C to be the measure with unit point mass at 0. Given U 
real numbers xi, X2, ■ ■ ■ , xjj and V non-real complex numbers, zi, Z2, ■ ■ ■ , zv, we define the 
measure r\ by 



d v(l) '■= ^2 a u dS(j -x u )+^2 b v ( dS (l - z v) + dS(j - z v )), 



1=1 



= 1 



where a\, az, . . . , au and &i, 62, ■ ■ ■ , W are indeterminants. It does no harm to assume that 
U and V are both greater than N. By reordering and renaming the x, z, a and b we will 
also write 

T 

t=i 



where we define 



d5{j-y) 



dS(rf-y) iiyeR; 
dSd - y) + dS(-y -y) if y e C». 



Clearly T = U + V. 

As wc alluded to previously, the proof of Theorem 4 relics on expanding Z w ^ x+ri ^ / Z w 
in two different ways and then equating the coefficients of certain products of ci, c 2 , . . . , Ct- 
One of the expansions of Z w ( x+71 ^ /Z w comes from Theorem 3, while the other comes directly 
from the definition of the partition function. 



Lemma 5. 



Pf (J + [^c^K N (y t , y t ,)]); t,t' = 1,2, ... ,T, 



(6.1) 



where J is defined to be the 2T x 2T matrix consisting of 2 x 2 blocks given by 



J := 



St.t' 



1 

-1 



t,t= 1,2,..., T. 



Remark. The Pfaffian which appears in (6.1) is an example of a Fredholm Pfaffian. This is 
the Pfaffian formulation of the notion of a Fredholm determinant and is discussed in [14]. 



A. Borodin, CD. Sinclair 



11 



We defer the proof of Lemma 5 and Lemma 6 until after the proof of Theorem 4. 
For each (I, m) and (L, M) we define the £, m, L, M-partial correlation function of 7 W to 
be Re, m ,L,M ■ R £ x C m -> [0,oo) where 

Re,m,L,M(x, z) : = 

y Q L ,M(xVa,zV (3) d\ L -t(a) d\ 2 (M- m ){P)- 

When i? == 772 == 0, we take x V a. — a and z V (3 = /3, so that Rq,o,l,m is a constant equal 
to ^l,m{Xl,m)- The partial correlation functions are related to the correlations functions 
by the formula 

i? £ , m (x,z) = — ^ Re,m,L,M(*,z), (6.2) 

(L.MJ 

and thus the partial correlation functions are one path to the correlation functions. Equa- 
tion (6.2) is still valid when £ = m = 0, though the constituent correlation 'functions' are 
actually constants. 

The partial correlation functions of 3>[ Gln ] of the special forms Rl,m,o,m and Ro,N,o,n 
have been studied by Akemann and Kanzieper in [1] and [12]. For general weight function 
w, the partial correlation functions of the form i?jv,o,n,o ar e (up to normalization) equal to 
the correlation functions of the (3 = 1 Hermitian ensemble with weight w. This connection 
will be exploited in Appendix A. 

Before stating the next lemma we need a bit of notation. Given non-negative inte- 
gers n and W, we define 3^ to be the set of increasing functions from {1, 2, ... ,77} into 
{1,2,..., W}. Clearly if W < n then 3^ is empty. Given a vector a e C w and an element 
u € 3^ , we define the vector a u € C™ by a u = {a u (i), a u (2), ■ • ■ , a u (n)}- 

Lemma 6. For each pair (L,M), 

n L . M (a, (3) d(X + rj) L {a) d(\ + n) 2M {(3) = 

L M , I m n 

\ Ii a uU)Y[K(k)\Re, m ,L,M(^n,^). 

<=0m=0„ e jl' oeJV < j=l fe=l > 

Remark. We will use the convention that 



This will allow us to keep from having to deal with the pathological correlation 'functions' 
Ro,o,l,m and i?o,o separately. 

Lemma 7. Suppose K : {1, 2, . . . , T} x {1, 2, . . . , T} -> R 2x2 is sucft tftat 

A-(i,0 = -^(t',t) T , 

and define K to 6e f/ie 2T x 2T antisymmetric matrix whose t, t' entry is given by 

K(t,t'). Then, 

T 

Pf[J + K] = l + ^ ^PfK t , (6.3) 



(L-^)!(M- 777)!2 M - 



L\M\2 M J RL Jc 
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where for each t £ dg, K t is the 2S x 25 antisymmetric matrix given by 

K t = [^(t( s ),t( s '))]; «y = i,2,...,5. 

Proof. This is a special case of the formula for the Pfaffian of the sum of two antisymmetric 
matrices. See [15] or [17] for a proof. □ 

Using these lemmas we may complete the proof of Theorem 4. First notice that Lemma 6 
and (6.2) imply that 



E E X I II^W) Tl b »(k) [^,m(x u ,Z B ). 
(l,m) ueiy 063^ 1 3 = 1 



k=l 



(6.4) 



This follows by summing both sides of the expression in Lemma 7 over all (L, M) and then 
reorganizing the sums over (L, M), t and m. 

Given u £ 3^ and <G we define (u: o) to be the element in 3^^ given by 



(uV 0)(n) := 



u(n) 
£ + o(n) 



if n < £ 
\in> I. 



Notice that each t e 3^ is equal to (u V o) for some u £ 3^ an d •= w ^ n £ + m = S. 
(This docs not preclude the possibility that either U or V equals 0, in which case t = D or 
t = u). It follows that we can rewrite (6.3) as 

pf(j+K)= E E E pfK (uv 0) . 

If we set K t j(t,t') = ^c t c t > K N (y t ,y t ,), then 



fe=i 



KN(x u (k), x u(k')) KN(x u (k), z o(n>)) 
(n) i **'u(fc / 



where fc and A;' are indices that run from 1 to £ and n and n' are indices which run from 1 
to m. Thus, Pf (J + K) equals 

t 



n^nw pf 



-K'Af(^u(fe) I 2 ; u(fc')) ^iv(a:u(fe),^ti(n')) 
(n) 5 *^u(fc' 



(6.5) 



and Theorem 4 follows by equating the coefficients of n ? =i a j ElfeLi ^fc m (6-4) and (6.5). 



7 Kernel Asymptotics for Ginibre's Ensemble 

We now turn to the large N asymptotics of the matrix kernel for Ginibre's ensemble of 
real matrices. In fact, we maintain our restriction to the case where TV = 2M is even, and 
consider the asymptotics of Kim as M — > oo. Throughout this section we will take <f> to be 
the function given by 



0(7) = cxp(- 7 2 /4 - 7 2 /4)^/crfc(v^|Im(7)|). 
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Notice that, since erfc(O) = 1, when 7 e R this reduces to exp(— 7 2 /2). 

The skew-orthogonal polynomials for this weight are reported in [9] to be 



2m 



7r 2m+ i(7)-7 2m+1 -2™ 7 : 



„2m-l 



7T2m(7) = 7 

with normalization {'K2m\^2m+i) v = 2v / 27r(2m)L The skew-orthogonality of these polyno- 
mials imply that 



S2M{l,i) 

DS 2 M{m') 



1 



2tt 
1 



M-l 

E 

m=0 
M-l 



2tt ^ 

m— 
M-l 



m=0 



7r 2m ( 7 )eA7r2m+l(7 / ) - 7r2m+l(7)gA7T2m(7 / ) 

(2m)! 

7T2m(7)^2m+l(70 ~ ^2m+l {l)^2 m (V) 

(2m)! 

eA7T2m(7)eA7f2m+l(7') - eA^2m+jX7^fA^2ro(^ 

(2m)! 



(7.1) 
(7.2) 
(7.3) 



The correlation functions are all of the form Pf K for an appropriate matrix K whose 
entries are given in terms of (7.1), (7.2) and (7.3). If D is a square matrix such that the 
product DKD T makes sense, then Pf (DKD T ) = Pf K • detD. And thus, if detD = 1 we 
have Pf K = Pf(DKD T ). That is, we may alter (and potentially simplify) the presentation 
of the Pfaffian representation of the correlation functions by modifying K in this manner 
by a matrix with determinant 1. When D is diagonal, the process of modifying K by D 
preserves the block structure of K. That is, the effect of modifying K by D affects changes 
at the kernel level and the correlation functions can be represented as the Pfaffian of a block 
matrix (c/. Theorem 4) with respect to a new matrix kernel Km dependent on Km and 
D. This will allow us to write the correlation functions of Ginibre's real ensemble in the 
simplest manner possible by 'factoring' unnecessary terms out of the kernel. 

It will be convenient to define cm, sm and eu to be the degree 2M — 2 Taylor polynomials 
for cosh, sinh and exp respectively. Explicitly, 



M-l 



C M (7) := Yl 



.,2m 



m=0 



(2m)! 



sm(7) == 



M-l 2m- 1 



; (2m- 1)! 



and 



cm(7) : = c m(7) + s m(7) 



2M-2 m 

^ m! 

m=0 



Theorem 7.1. The £,m- cor relation function of Ginibre's real ensemble of 2M x 2M ma- 
trices is given by 



i?,, m (x, Z ) =Pf 



K 2 M{Xj,Xf) K 2 M(Xj,Zk<) 
K 2 M{zk,Xj') K 2 M{zk,Zk') 



3,f -- 



1,2, 



J:. 



1,2,.. .,m, 



where 



#2m(7,7 ; ) = 

is given as follows. Let x, x' € WL and z, z' G C*, then: 



DS2M{l,l') 5 , 2 m(7,7') 
-S2m(i',i) IS 2 m{i,i') + £(7,7') 
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1. The entries in the real/real kernel, K 2 m(x,x'), are given by 

• S 2 m{x,x') = -j= — e xx euixx) H -j^^ru (x, x'), where 

\Z2ir \/2ir 



oM-3/2 / i 12- 

Tm{z,x') = J^y^z^ -l[M- -,- ) ; (7.4) 

• DS 2M {x,x') = = — (x' - x)e xx e M (xx'); 

\/2tt 



IS 2M {x,x') 

p -x 2 /2 f x' 2 /2 t _ p -x' 2 /2 f x 2 /2 t 



2^ 



■/2 rx I A e -t e -x "/2 rx ,z £ - t 

— sgn(x') / — pCM(xV2£) dt - — -= sgn(x) / — F CM(x / V2t)dt. 
7T Jo V* 2 V 7r Jo Vt 



2. The entries in the complex/ complex kernel, K 2 m{ z , z '), ar e given by 
-a (*--') 

/27T 

e 2 



S 2M {z,z') = 16 2( * Z) (z 1 - z) J eric (V2Im(z)) crfc (V2lm{z'))e- zz ' e M {zz' 



• DS 2M {z,z') = = — (z' - z)J eric (V2Im(z)) eric (V2rm(z'))e- zz 'eM(zz')] 



_ _ -\{-~-') 2 i _, 

IS 2M {z,z') = — {z'-z) Verfc(v / 2Im(z))erfc(v / 2Ini(z'))e- zz e M (^')- 

a / sir * 



/2tt 

5. T/ie entries in the real/complex kernel, K 2 m{x,z), are given by 
• S 2 m(x,z) = = — (z- x)y eric (V2lra(z))e~ xz eM(xz); 



_ g— 5(2: — z) 2 / 

5 2 m(z,2:) = -j= — Jeric(V2Im(z))e~ xz e M (xz) 



+- 



2ir 



\J eric (V2lra(z))r M (z,x); 



• DS 2M (x,z) = = — (z - x) J eric (V2Im(z))e xz e M (xz); 

y2ir 

. e -\(x--,) 2 I — 

• IS 2M {x,z) = -= — y crfc (V2Im(z))e~ xz e M (xz) 



e * 
+- 



■-y/erfc (V21m(z))rM(z,x). 



2ir 

Theorem 7.1 allows us to derive the M — > 00 limit of K 2 m- 

Corollary 7.2. Let x and x' be real numbers, and suppose z and z' are complex numbers 
in the open upper half plane. We define K = lim K 2 m- Then, 

M^oo 
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1. The limiting real/real kernel, K(x,x'), is given by 



K(x,x') 



L {x' - x)e-^ x - x '? 



J_ -i(x-x') 2 
/2tt 



| sgn(a; - a;') erfc (^=M) 



TTie limiting complex / complex kernel, K(z,z'), is given by 
K(z,z r ) = J erfc (V2Im(z)) erfc (V2Im(z)') 



' (z' - z)e-^ z - z '^ i(z' - z)e-^ z -~"> 2 
i(z' - z)e-^ T '- z 'f -{z 1 - z)e-^-----'f 



3. The limiting real/complex kernel, K(x,z), is given by 
K(x,z) = -±=yJcrfc(V21m(z)) 



[z — x)e ii x z ) 2 i(z — x)e *( x z ^ 
g— -^{x— z) — ie~ 2( x ~ z ) 



We now turn to the kernel of Ginibre's Real Ensemble in the bulk. The circular law for 
N x N matrices with iid Gaussian entries says that, when normalized by N^ 1 / 2 the density 
of eigenvalues becomes uniform on the unit disk as N — > oo (See [11] for a proof of this fact 
when the entries are iid Gaussian, and [2] and [18] for more general results). This gives us 
the appropriate scaling when considering the matrix kernel in the bulk. Specifically, in this 
section we will be interested in the large M limit of K 2 m(uV2M + s, u\j2M + s') where u 
is a point in the open unit disk, and s and s' are complex numbers. 

When it is real we expect that the limiting kernel under this scaling should yield K(s, s'); 
indeed this is the case. When u is nonreal a new kernel arises. 

Theorem 7.3. Let — 1 < u < 1 be a real number, let n, r 2 , ■ ■ ■ ,re € R and S\, s 2 , ■ ■ ■ ,s m be 
in the open upper half plane. Set, 

Xj — u\J 2M + rj j = l,2, ...,£; and z k = u\Fl~M + Sk fc=l,2, ...,m. 

Then, 



lim R( m (x,z) = Pf 

M^oo 



K(rj,rj>) K(rj,s k >) 
K(s k ,r r ) K(s k ,Sk<) 



j,f = 1,2,...,£; 
fc, k' = 1, 2, . . . , m, 



where K is given as in Corollary 7.2. 



Theorem 7.4. Let u be in the open upper half plane such that \u\ < 1 and suppose 
si, S2, • • • , s m e C. Set, 



Zk = u\ 2M + Sk k = 1, 2, . . . ,ra. 



Then, 



lim Rq to (— , z) = det 

M— >oo 



1 / l*fcl 

-exp -- 

7T 



\Sk> 



+ S k S k i 



k,k' = l 
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Remark. The function 

( >\ 1 f |s|2 |s '' 2 + ~\ 

K (s, s ) — — exp I — + ss' I 

is the limiting kernel of Ginibre's complex ensemble. Thus, the limiting correlation functions 
in the bulk of Ginibre's real ensemble off the real line is identical to the limiting correlation 
functions in the bulk of Ginibre's complex ensemble. 



8 Proofs 

8.1 The Proof of Lemma 5 and Lemma A. 3 

In the case of the real asymmetric ensembles Y — C and b = 1. In the case of the Hermitian 
ensembles Y = K and b = y/]3. We start with 

(qn\q n ') w (\+ri) = (q n \q n ')\+7, = {qn{i)z\+riq n '{i) - £\+riq n {i)qn'{i)) d(\ + n)(-f). 

An easy calculation reveals that this is equal to 



(q n \q n ')wx + -r y^ct(qn(yt)e\qn> (yt) - qn'(yt)e\q n {y u )) 



t=i 



T T 



_ ^ c t c tl q n (y t )q„>{yt>) £(y t ,yt>)- 

t=i t'=i 

Next we define A to be the Nb x 2T matrix given by 



A := 




n = 0,l,...,Nb-l and t=l,2,...T, 



and the 2T x 2T matrix B, 
B := — J + 



y/c t c t > £,(y t ,yt>) 




t,t' = 1,2,...,T. 



We define C to be the Nb x Nh matrix given by C = (Uq ) ; the n, n' entry of C is n n ,n'- 
A bit of matrix algebra reveals that 



Z w(\+v) _ P f (c-t - ABA T ) _ Pf(B- T - A T BA) 



Z w 



Pf(C- T ) 



Pf(B- T ) 



(8.1) 



where the second equality comes from the Pfafhan Cauchy-Binct formula (sec Appendix B). 
Notice that 



B T 



i,i' = l,2,...,T. 



J- 

[0 y/ctCf E(y t ,y t ')_ 

And from Lemma 7, Pf(B) T = Pf(— J) = (— 1) T . A bit more matrix algebra reveals, 



A T CA = 



y/c t c t > DS N (y tl y t >) ^c t c t > S N (y t ,y t >) 
-y/c t ct> S N (y t ',y t ) y/ctc t/ IS N (y t ,y t >) 



t,t' = 1,2,... T. 
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Using these facts and simplifying (8.1) we find 

zw = (-if Pf ( - J - WwK N (y t ,y t ,)}) ; t,t' = 1,2, ... ,T, 

and the Lemma follows by using the fact that if E is an antisymmetric 2T x 2T matrix, 
then Pf(-E) = (-1) T Pf(E). 

8.2 The Proof of Lemma 6 and Lemma A. 4 

We start with 

l!j J, 2 m J J M VlMccP) d(X + r,) L (a) d(\ + V ) 2M (f3). (8.2) 

Notice that in the case of the Hermitian ensembles that this is equal to Z w ^ +V ^ when L = N 
and M = 0, and the proof of Lemma A. 4 follows from the proof recorded here by setting 
every instance of M to 0. 
First we write 

L L 

d(X + v ) L (a) = Yl (rfAi(a,) + d m ( aj )) = E E dve(a u )d\ L -e{cx u >), (8.3) 

where given t € 3^ , we define t' to be the unique element in 3^_ n whose range is disjoint 
from t. Notice that since u' appears in the summand on the right hand side of (8.3), the 
inner sum is not actually empty when I = 0; in this situation the summand is equal to 
d\ L {a). 
Similarly, 

M 

d{\ + rj) 2M (P) = E E ^2m(/3 B )rfA 2 ( M -m)(/3 ')- 

"1=0 ve3™ 

Thus, (8.2) equals 

EE E E L\M\2 M / / Q L,M(a,l3)dri t (a u )d\ L - e (a u ,)dri2 m (P )d\2(M- m )(Pt,>)- 



We can relabel the a and (3 in the integrand in any manner we wish, and in particular we 
may make the integrand independent of u and 0. In particular, if we set i e 3™ to be the 
identity function on {1, 2, . . . , n}, and since the cardinality of 3™ is (^), we find that (8.2) 
is equal to 

L M 

EE 

£\(L — t)\m\(M — m)\2 M 
x J J {J Re J c n L M( a ,f 3 ) d Ve(ai)di 1 2, n {f3 i )^d\ L ^i(a i ,)d\2(M- m )(P l >)- (8-4) 

R L-t C M- m 

Now, 

l £ U 

dr]e(cxi) = ~[[vi( a j) = nE "^^' ~ x u)- 
j=i j=i tt=i 
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We may exchange the sum and the integral on the right hand side of this expressions by 
using the set, ^ of all functions from {1,2,..., £} into {1,2,..., U}. Specifically, 



f 1 



and similarly, 



oe$¥. ^ fe=i J 



t>e5. 

Thus, 

/ / n LM (oc,f3)drji(a i )dr]2 m (f3 i ) 
r £ m ) 

= Y Y \ Il a »U) n 6 "( fc ) [2 m O iiM (x u Va iS z V/3 t ,). (8.5) 

Notice that if u or D is not one-to-one then | A(x u V oh, z„ V /3 t ) | = 0. We may consequently 
replace the sums over ^ and ^ with their respective subsets of one-to-one functions. 
Moreover, since SIl,m is symmetric in the coordinates of each of its arguments, we may 
replace each one-to-one function in these sums with the increasing function with the same 
range so long as we compensate by multiplying by l\ and ml. Lemma 6 follows from the 
definition of Re, m ,L,M by substituting (8.5) into (8.4). Lemma A. 4 follows from the fact 
that R n = R n ,o,N,o/Z w . 

8.3 The Proof of Theorem 7.1 and Corollary 7.2 

It shall be convenient to introduce the following variants of S2M and DS2M' 

£ , n 7r 2m(7) e A7f2m+l(7') - 7r 2m+l(7) e A7T2m(7') 

S 2 m(7,7 ) ■■= L (a^jj ' 

m=0 v ; 

and 

"HO / l\ . 7T2m (7)7r2 m +l(70 ~ ^+1(7)^(70 

^2M(7,7 ) •= ^ ^ , 

m=0 y ' 

The following lemma gives a closed form for these functions. 

Lemma 8.1. Let x be a real number, and suppose z and z' are complex numbers. 

1. S2m{z,x) — <p(x)eM{zx) + rM(z,x), where ru is given in (7-4); 

2. DS 2M (z, z') = (z 1 - z)e M {zz'). 
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Proof. First we compute e\ir 2m and cx^m+i- We start by noticing, 

1 f°° 1 f x 1 f°° 

^g( x ) = 2 y ) s § n (y - x ) d v = -2 J 9{y) dy+ 2 J 9 ^ dy ' 

When g(y) = e~ y2 / 2 y n , we may evaluate the latter two integrals in terms of the incomplete 
gamma functions. 



n(n—l)/2yf Tt+1 x 2 
L 1 ^ 2 ' 2 

We immediately conclude that 



_ 2 ("-i)/2 sgri ( x ) . 7 (n±I, if n is even; 

if n is odd. 



e\W2 m (x) = -2 m - 1/2 sgn(a;) • 7 + 1, 



and 



eA7T2m+lO) = 2 r ' 



r m + 1,- 



-mrf 



(8.6) 



(8.7) 



where in the second equality we used the fact that T(a + l,x) = aT(a, x) + x a e 
Using (8.6) and (8.7), we may write 



S 2 m{z,x) = cj>{x)c M {zx) 

(M-l 

+ s g n 0) I 



M ~ 1 2 m ~ 1 / 2 



-J (2m)! 

m=0 v ; 



2\ M -i 2" 1 - 1 / 2 



2' 2 7 ^ (2m- 1)! 

' m=l x ' 



Next, we use the fact that 

7(a + 1, .t) = 07(0, x) — x a e~ x , 
so that the second sum in this expression becomes 



M-l 

E 



2(m-l)-l/2 

^(2(m-l))f 
Consequently, 



,2(m- 



i)+i 7((™-i)+^)+^x)5: 



2' 2 



SW(z,a:) = <j>(x)(c M (zx) + sgn(x)s M (z\x\)) + 



2M-3/2 
(2M - 2) 



= ^(i)e M (2i) + 



2M-3/2 
(2M- 2)! 



sgn^^SlM-^, 



^ (2m -1)! 

1— 1 v y 



Sgn^z^-^IM--,^- 



1 x v 



1 x 2 
2' 2 



2' 2 
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Turning to DS 2 m, 



z 2rn( z ,2rn+l _ 2mz ,2 m -^_( z 2rn+l _ 2mz 2rn-ly2rn 

DS 2M (z, z ) = ^ — ~- { 

m=0 v ; 

M_1 .2m,l2m+l 7 2m+l J2m M_1 2ro-l J2ra 7 2mJ2m-l 

Ez z z z ^ — ^ z z z z 

n (2m)! + ^ (2m- 1)! 

m=0 v ' m=l v 7 

{ M - X z 2m z l2m ^m-l^m-l ~| 

= (z' — z)e M (z'z). □ 

With a closed form for SW and DS 2 m in hand, we are ready to prove Theorem 7.1. 
Proof of Theorem 7.1. From Lemma 8.1 we have that 

r, , I \ <t>{x) ^ , l\ 4>{x)4>{x') , (j)(x) ,, 

S 2 M(x,xf) = ^=LS 2M {x,x) = yv ; 11 '- e M {xa!) + -j=Lr M {x,x) 

V27T V27T V27T 



and 

D^2m(a;, a; ) = -= — DS 2M (x, x') = ' - (a; - x)e M (a;x ). 

V27T V27T 



Now, 

0(x)0(x') = e^O^' 2 ) = e -|(x-x') 2 e -xx' ; 



g 2 X ) / (^>f 

S 2M {x,x') = -= — e~ xx e M {xx') + -±=Lr M {x,x') 



and therefore, 



and 



The computation of IS 2 m(x,x') is a bit more involved. From (7.3), (8.6) and (8.7), we 

see 

M-l nm , , , 2 - 



DS 2M (x,x') = = — (x' - x)e xx 'e M (xx') 

V27T 



v m=0 v ' v ' 

1 ^ 2 m / 1 x 2 \ . . I2m _ x ,2 /2 
"2V5fE (2^)!-H m+ 2'Tj Sgn(x)X 6 



IS 2 m{x,x') is clearly skew-symmetric in its arguments; looking at the first sum in this 
expression we thus find, 

M-i om „ r x' 2 /2 



V -=— -sgn{x')x 2m e- x / 2 / r-We-tdt 
2^F^ (2m)! 7 

1 



/■x' 2 /2 e -t (M- 1 9 m 1 
^e"* /2 sgn(x') / e — { Y -j^-:t m x 2m dt, 
A 'Jo sft \„t- (2m)! f 
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where on the left hand side we have replaced the lower incomplete gamma function with its 
integral definition, and on the right hand side we have exploited the linearity of the integral. 
The sum on the right hand side of this equation is equal to c M (xy / 2~i) 7 and thus 



IS 2M {x,x') = 

sgn(a/) / —= c M {xV2t) dt sgn(x) / —=c M {x'V2t)dt. 

Jo yt W 71 Jo yt 



Turning to the complex/complex entries of K 2 m, if z is assumed to be in the open upper 
half plane then e\Tr n (z) = i-K n ^z). From this we see that S 2 m(z, z') = i(f)(z)<j)(z')DS 2 M(z, z'), 
DS 2M (z,z') = <j>(z)<j>{z')DS 2M {z,z') and IS 2M (z,z') = -<j>(z)<f>(z')DS 2M (z,z'). 

Next, we define 

iP(z) = ei<- z2 -~ 2 \ 

Notice that 
and thus 

4>(z)4>{z') = iP(z)ii(z')e-^ z - z ' )2 ^erfc (V2 Im(z)) erfc (V2 lm(z)')e~ zz ' . 
Using this and Lemma 8.1, we conclude that 

S 2M (z,z') = 4>(z)4>(z') te ' {z - z)yerfc (y/ZSm(z)) erfc (V2lm(zY)e~ zz e M (zz'), 



2tt 

DS 2M (z,z') = i)(z)i)(z') e 2< 1_ Z) (z' - z)Jerfc (V2\m(z)) erfc (V21m(z)')e-"' e M (zz'), 

\/2ir 



and 

_ -i(z-z') 2 _ i _ 

S 2M (z,z') = i){z)il>{z') = {z' - z) J erfc (s/2lm(z)) erfc {V2lra{z)')e~ zz e M {zz')- 

\/2ir 

Lastly we look at the real/complex entries of K 2 m- As in all other cases, 
DS 2M (x,z) = M^DS 2M (x,z), 

and it is easily verified that 



4>(x)4>(z) = ^(z)e'i^ z ^yJer{c(V2lm(z))e- 



Thus, 

e -i(x-z) 2 1 — 

DS 2M {x,z) = ip(z) — — (z - x)ycrfc (V2 lm(z))e~ xz e M {xz). 



/2tt 
Since, 

a 1 \ .<t>( x )<t>( z ) -jr^ </>(z) ? , > 

S 2M {x, z) = 1 — DS 2M {x, z), S 2M {z, x) = —=S 2M \z, x). 

v2tt v2tt 
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and 



It follows that 



IS 2 m{x,z) = -i^=S 2 M{z,x), 



S2m(x, z) = ip(z) == — (z — x) J eric (V2Im(z))e~ xz eM{xz) 1 



2tt 



i(z-z) 2 



2tt 



erfc (v / 2Im(z))e xz e M (a;z) 



+ - 



e * 



2tt 



y^erfc (v / 2Im(z))r M (2, x) 



and 



e -5(^) 2 / — -__ 

IS 2M (x,z) = ip(z) \ j= — A/crfc (V2Im(z))e xz e M (xz) 



+ - 



erfc (v / 2Ini(z))r M (z, x) 



Clearly, ip{x) = ip{x') = 1, and thus we find that 

K N {i,i) 



Mi) o 

V(7) 



It follows that, if we define K to be the matrix 
K = 



\p(i) 

V(7') 



K N (xj,Xf) K N {xj,z k >) 
K N (z k ,Xj>) K N (zk,Zk>) 



j,f = 1,2,... A 

k, k' = 1, 2, . . . , m, 



and D to be the diagonal matrix 

D = diag(V>(zi),^(Sl), . . . , ip(xi),ip(xj),ip(zi), V>(zT) • • • , V'C^m), V'C^O), 

then 

i?,, m (x,z) =Pf(DKD). 
But, since ^("z) = VK- 2 ) -1 , we have that detD = 1, and i?^ m (x, z) = Pf K as claimed. □ 
Proof o/ Corollary 7.2. We first make use of the fact that 

lim ejf(z) = e z 

pointwise on C. This simplifies all terms in the kernel except the IS 2 m term. It remains to 
show that rM(z,a;)-tOasM^oo, and that 



\ sgn(x - x 1 ) + lim IS 2M {x, x 1 ) — \ sgn(x - x') erfc ( — ^ 



(8.9) 
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The first of these facts is easily seen by noting that j(M - 1/2, x 2 /2) < T(M - 1/2), and 
hence 



\r M (z,x)\ < 



\ Z \2M-1 2 M-3/2 r(M _ 1/2) \ Z \2M-1 2 ~ M + 3 / 2 ^ T(2M - 2) 



T(2M - 1) 



T(2M-1)T(M- 1) 



where we arrive at the latter equality we have used Legendre's duplication formula on the 
T(M- 1/2) term. Thus, 



\r M {z,x)\ < 



|2M-1 2 -M+3/2 ^ 



r(M-l) (2M-2) 



and it is easy to see that both terms in this product go to as M — > oo, independent of the 
value of z. 

To establish (8.9) we start with 



I M {x,x') :-- 



Vt 



CM{x\2t) dt 



5.10) 



I(x,x') := lim Im(x,x') 

M^oo 



Since the terms in cm are all positive, from the Monotone Convergence Theorem, 

■ cosh(xv / 27) dt 

Jo Vt 

/ It' I + t\ 

erf 



2 



The latter equality follows from the fact that 



V2 



erf 



\x'\ -x 



V2 



(8.11) 



o y/t 



cosh(aVi) = ^ e » 2 /4 J erf g + ^ _ erf _ ^ 



which can be verified via differentiation. 
Now, 

lim IS 2 m{x, x 1 ) = —^-= le~ x2/2 sgn(x')I(x, x') - e~ x ' 2/2 sga(x)I(x' , x)\ 

M— >oo 2\JTT I J 

1 



sgn(x') erf 



\x'\+x 



V2 

— sgn(ai) erf 



sgn(x') erf 
x' + \x\ 



V2 



|x'| -x 
+ sgn(x) erf 



x - m 

V2 



= --erf 



V2 



It follows that (8.9) can be written as 

1 „/x-x'\ 1 



- sgn(x - x') - - erf 



V2 



■ sgn(x — x 1 ) — — sgn(x — x') erf 



\x-x'\ 

~7T~ 



where we have exploited the fact that erf is an odd function. We arrive at the form for 
IS2M stated in the corollary using the fact that crfc = 1 — erf. □ 
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8.4 The Proof of Theorem 7.3 and Theorem 7.4 

In order to prove Theorems 7.3 and 7.4, it is necessary to investigate the asymptotics of the 
partial sums the exponential function. 

Lemma 8.2. Let u^±l be a complex number, and let (v m )^ =1 be a sequence of complex 
numbers satisfying 

VM = u 2 + O (i/Vm) as M — > oo. 

Then, as M — ► oo, 

p -2(l-u 2 ) p 2M(l-u 2 )„4M 

e- 2M ™e M (2Mv M ) ~ 1 - ^ ■ - 7=^—. 

2ttu 2 (1-u 2 ) s/m 

In particular, when — 1 < u < 1 is real, 

lim e- 2Mt,M e M (2Mv M ) = 1. 

M^oo 

Proof. Set u = u M . We start by writing 2¥« = v(2M - 2) + 2t>. Thus, 



-2Mb 

e e M 



(2Mv) = cxp (-(2M - 2) + jj^y)) eM ((2M - 2) + 



We write 

Clearly w M = m 2 + 0(l/\/M). Under this hypothesis, and since m 2 7^ 1 , equations (2.9), 
(2.15) and (1.7) of [4] imply that 

p -2(l-tt 2 ) „2M(1-m 2 )„4M 

e -m-2 )w {{2M _ 2)w) „ 1 _J_ . £ n 

2iru 2 (l - u 2 ) y/M 

The second statement of the lemma follows from the fact that if u is real and \u\ < 1, then 

e 2M(l-« 2 ) w 4M _ e 2M(l-u 2 +21og|u|) 

and 1 — u 2 + 2 log |u| is negative when u is in (—1, 1). 

Lemma 8.3. Suppose u,sel with \u\ < 1 and suppose s' e C. Let 

oM-3/2 

^ , /rr~r~r ;n2M-1 



i? M = exp (-Mm 2 - uV2MRe(s')) (2M _ 2) | ("^ + s ') 



(2M - 2) 

x 7 f M - Mm 2 + uV2Ms + y ] • 



Tften, lim J?m = 0. 

M— >oo 

Proof. Clearly, 



7 (m - 1, m„* + uvms + £) < r (m - i) . ^ 2 -»««Ee^2> , 
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where the second equation follows from Legendre's duplication formula for T. Now, 

luV^M + s'l™- 1 ~ | M |2M-l 2 M-l/2 M M-l/2^ 

Thus, 

\Rm\ 



M M+l/2 

^| u |2M-i ^_ exp (_ Mu 2^ exp (- _ uy /2M 8 y 



M\ 

Using Stirling's approximation for Ml, we find 

/2 

\Rm \ ~ ^-|u| 2M_1 exp (M(l-u 2 )) exp(-it\/2Ms) 



V2 
2\u\ 



exp (M(21og|u| + 1 - u 2 )) exp ( - mv / 2Ms). 



We observe that 2 log |tt| + 1 — u 2 < since |u| < 1, and hence R(M) — > as M — > oo. □ 

Proof of Theorem 7.4- Let u be a point in the open upper half plane with modulus less 
than 1, and suppose s and s' are complex numbers. For all but finitely many values of M, 
z = u\J1M + s and z' — u\j2M + s' arc in C*. Thus, in this case, we need only consider 
the asymptotics of the complex/complex kernel under these substitutions. 
We will make use of the fact that if a: is a real number, 



erfefa;) = — =r f -,x 



y/ir\x\' 



Consequently, 



—p " — 1= — exp ( - 2MIm(u) 2 - 2v / 2MIm(u)Im(s) - Im(s) 2 ) , 

erfc v/2Im(W2M + s) ^ J > \ > LU. ( 8 .12) 

^/2Im(u) <J M-k 



Now, by Theorem 7.1, 

^2m(mV21?+s, U \/2M + s') = ^ ~jh -i( s - s '^ 

V 27T 



x yj erfc (V21m(uV2M + s)) erfc (V2Im(u\/2M + 

x cxp(-2Mu 2 - ?iV2M(s + s') - ss') e M (2Mu 2 + uV2M(s + s') + ss') . 



Therefore, by Lemma 8.2 and (8.12), 

DS 2M (uV2M + s, uV2M + s') ~ ifl^e-^-*') 2 



2tt 

e -4MIm(«) 2 e -2\/2MIm(u)(Im(s)+Im(s')) e -Im(s) 2 -Im(s') 2 



X 1- 



-2(l-tt 2 ) 



2Im(u)VM7r 

e 2M(l-tt 2 ) u 4M~ 



2ttm 2 (1 -u 2 ) 



M 
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It is easily seen that 



-4MIm(u) 2 -2\/2MIm(ti)(Im(s)+Im(s')) 

lim — = 0, 

M->oo 2Im(u)V 'Mir 



and, 



e -4MIm(u) 2 e 2M(l-« 2 ) u 4M 



2M(l-|u| 2 +21og|u|) 



Since \u\ < 1, we have 1 — \u\ 2 + 2 log |u| < 0, and therefore 



lim e -2V^MIm( U )(Im( a )+Im( S ')) 



e -4MIm(u) 2 e 2M(l-u 2 ) u 4M 



= 0. 



We conclude that 



lim DS 2M (uV2M + s, u\p2M + s') = 0. 

M^oo 



And, since IS2m(z,z') — —DS 2 m{z, z'), 



lim IS 2 m(uV2M + s, uV2M + s') = 0. 



Turning to S 2 m(uV~2M + s, u\/2M + s'), we set 

Vm{s) = exp (-2iv / 2MIm(u)Re(s)) . 
From Theorem 7.1, we have 



t] M (s)vm(-s')S2m{u\/2M + s, uV2M + s') 



2iV2Mlm(u) + (s - s') 



x exp (2iV2MIm(u)Re(s)) exp ( - 2iv / 2MIm(u)Re(s / )) 
x exp ^4M(Im(w)) 2 - 2iV2MIm(u)(s - s') - l(s - s'f^J 

x ^erfc |2\/Mim(w) + \/2Im(s)| erfc |2-/MIm(u) + V2Im(s')| 

x exp {-2M\u\ 2 - {su + s'u)V2M - ss'^j e M (2M\u\ 2 + (su + s'u)V2M + ss'^j 

From Lemma 8.2 and (8.12) we find that 



T] M (s)r]M(-s')S2M(uV2M + s, uV2M + s') exp --(s - s') 2 - Im(s) 2 - Im(s') 



7T 
1 

— exp 

7T 



\s\ 2 W\ 2 



2 2 

Next, we set D to be the 2m x 2m diagonal matrix given by 

Dm = diag(?7M(si),f7M(-si), . . . ,r)M(s m ),r)(~s M )) 



.S'.S 
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noting that detDAf = 1. It follows that 
Jim i? .m(-,z) = lim Pf (D M [K 2M (s k , s k >)] 7 ji .,_ 1 D^ f 1 

I exp (- 

1 



= Pf 







J k,k' = l 



= dct 



■ exp 



Isfel 2 |s fc /| 2 



k.k' = l 



where the last equation follows from Section 4.6 of [16]. 



□ 



Proof of Theorem 7.3. First we consider the case where s and a' are both in the open upper 
half plane. From Theorem 7.1, 



S 2M (2M)(uV2M + s,uV2M + s r ) = 



i e -5(«-»') 



(a' - s) J eric (\/2Im(s)) erfc (V2Im(s')) 



x exp ( - 2Mu 2 - uV2M(s+f) - ss')e M (2Mu 2 + uV2M(s + s') + ss'), 

, , e -i(s-s') 2 i 

DS 2M (2M)(uV2M + s,uV2M + a') = j== — (s' - sWerfc (\/2Im(s)) erfc (v^Im(s')) 



27T 



exp ( - 2Mu 2 - uV2M (s + s) - ss') e M (2Mu 2 + uV2M (s + a') + ss') 



and 



IS 2M {2M)(uV2M+8,uV2M+a') = 



_ e -Jr(s-s') 2 



2tt 



(a' - a) ^crfc (V2Im(s)) erfc (V2Im(s')) 



x exp ( - 2Mu 2 - uV2M (a + s') - ss')e M (2M« 2 + uV2M (a + s') + ss') . 

By Lemma 8.2, these converge to the appropriate entries of the complex/complex kernel 
K(s,s') as M —> oo. 

Next we turn to the case where u and r are real, and s is in the open upper half plane. 
In this case, Theorem 7.1 yields, 

S 2M {uV2M + r, uV2M + s) = - 



2tt 



(s-r)^/erfc(V2Im(s)) 



x exp ( - 2Mu 2 - uV2M(r + s) - rs)e M (2Mu 2 + u\p2M{r + s) + rs), 



S 2M (uV2M + s, uV2M + r) 



e 2 



2tt 



erfc (V2Im(s)) 



+ 



x exp ( - 2Mm 2 - uV2M(r + a) - rs) e M (2Mu 2 + u\p2M{r + s) +rs), 
^L^erfc (V2Im(s))e-3(" 2 +« 2 ) exp ( - Mu 2 - uV2MRc(s)) 



2M-3/2 



sgn(wV2M + r){uV2M + s) 2 * 1 "^ M - -, Mu 2 + «V2Mr 



(2M- 1)! 
and thus, by Lemmas 8.2 and 8.3, 

~ ie~ i(' — ^ 2 / 

lim S 2M (uV2M + r, uV2M + s) = = — {s-r)J erfc (\/2Im(s)) , 

M— >CX) V27T V 
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and 



_ 1 (r-s) 2 

lim S 2 m(uV2M + s, uV2M + r) = 6 ' 



^/erfc(v / 2Im(s)). 



The limiting values for DS 2M (uV2M + r,uV2M + s) and IS 2M {uV2M + r,u\/2M + s) 
follow from this as well, since 

DS 2M (z,z') = -iS 2M {z,~z) and IS 2M (z,z') = iS 2M (z',z). 

Finally, we turn to the case where u, r and r' are all real. Here, Theorem 7.1 implies 



that 

S 2M {uV2M + r, uV2M + r 1 ) 



2ir 



x exp ( - 2Mu 2 - uVIM (r + r 1 )- rr') e M (2Mu 2 + u\Fl~M (r + r') + rr') 

I 1 2 

+ —=e~z r exp ( - Mu 2 - uV2Mr) 

2M-3/2 



(2M - 1) 



— - sgn(uV2M + r'){uV2M + r) 2M "S ( M - ^Mu 2 + uVWlr' + — J . 



and 



DS 2M (uV2M + r, uV2M + r') = ^ _^ e-^-r'f 

\/2ll 



x exp ( - 2Mu 2 - uV2M(r + r') - rr')e M {2Mu 2 + uV2M(r + r') +rr'), 
From Lemmas 8.2 and 8.3, we see that 

lim S 2M (uV2M + r,uV2M + r') = — Le^-'") 2 , 

M-OO y/^ 

and 

lim DS 2M (uV2M + r,uV2M + r) = ^ ZlK ~^ r - r "> 2 . 



\r-r'\ 
V 



All that remains to show is 

lim IS 2M {uV2M + r, uVIM + r') = \ sgn(r' - r) erfc 

M^oo 2 

First we write 



IS 2M (x, x 1 ) = sgn(x') ^ I(x, x') - / -^C M (xV2t) 

2 V 7T 1 Jo Vt 



dt 



2^ 



■sgn(a;) \ I(x',x) 



f 

Jo 



2 1 2 e -* 



V 



C M {x'V2i)dt 



where I(x, x') is given as in (8.10) and Cm = cosh —cm- That is, 

~ er x2/2 e~ x ' 2/2 
IS 2 m(x, x 1 ) = -j= sgn(x')/(x, x 1 ) _ sgn(x)I(x' , x) 



+ 



2^ 

■ sgn(a;) 



2^ 



2 1 2 e~* 



-x 2 /2 



r C M {x'Vlt)dt 



sgn(x') 



*' 2 ' 2 e-* 



Vt 



C M {xV2t)dt. 

(8.13) 
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Making the substitutions x = u\J1M + r and x' = uv / 2M + r', and assuming that M is 
sufficiently large, (8.11) yields 



J im 7= sgn{x')I{x,x') 



icrfc(^) ifu>0; 



Consequently 



,l im \ o 7= sgn(x')I(x, x') - sgn(x)/(x', x) 



— erfc [ — =— | — - erfc 



V2 / 4 V V2 



r — r 



1 c (r — r 1 
■— crt 



2 V V2 

Thus, if we can show that the second line of (8.13) goes to as M — > oo, we will have 

fl 1 1 \ ( r - r' 

m^o\ 2 sgn ( x ~ ^ ) + /5< 2M(a;,a;')| = g s S n (r - /) - - erf I — 

= - sgn(r — r ) eric ' 



as desired. 

We thus consider 

e - x ' 2 / 2 r x ' A/2 



2V^ 

Clearly, for any v > 0, 



— sgn(a;) / -^C M {x' 'V2t) dt. 

* JO Vt 



ro=2M-l v y 



Thus, 

r x 2 /2 t _ r x 2 /2 t 



/ ^C M (x'V2l) eft < / — = V^erf I i-L ) < VtF. 
Jo \t Jo vt Vv2/ 



It follows that 



Jo \t 



Jo Vt ~ 2 ' 

and making the substitution a;' = uv"2M + r', this goes to as M — > oo. □ 

Appendix A Correlation Functions for (3=1 and (3 = 4 
Hermitian Ensembles 



In this appendix we will use the Pfaffian Cauchy-Binct Formula (see Appendix B) in order 
to derive the correlation functions of the (3—1 and (3 = 4 Hermitian ensembles. We will 
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keep the exposition brief, but will introduce all notation necessary for this appendix to be 
read independently from the main body of the paper. We reuse much of the notation from 
main body of the paper so that we may also reuse the same proofs. For convenience, N will 
be a fixed even integer; similar results are true for odd integers. 

Given a Borel measure v on R we define the associated partition function to be 



Z V :=^J |A( 7 )|d^( 7 ), 

^ ■ JR« 



where A (7) is the Vandermondc determinant in the variables 71,72, ■ ■ • ,7at and vn is the 
product measure of v on R N . When (3 = 1 we define the function £ : R 2 — > {— ^, 0, 5} and 
the operator e„ on L 2 (v) by 

8(7,7') := 77Sgn(7-7') and £^(7):=/ g{y)£{y,j)dv(y). 

When [3 — 4 we define £(7,7') := and e l/ g(y) := g'(y). We use e„ to define the skew- 
symmetric bilinear form on L 2 {v) given by 

(g\h)v~ / (g(-f)e„h(j) - £v9(j)Hj))dv(j). 
Jr 

Theorem A.l. Let b := \f]5, and let q be a family of Nb monic polynomials such that 
deg<7„ = n. Then, 

where = [{q n \q n >)v\; n,n' = 0,1, ... , Nb — 1. 

This theorem follows from de Bruijn's identities [6]. 

We set A to be Lebesgue measure on R. If there is some Borel measurable function 
w : R — > [0, 00) so that v = wX (that is, dv/dX — w) then we define Z w := Z v . Clearly, 

Z w = ^J^n N ( 1 )dX N ( 1 ) where ^(7) := | J] w{ ln ) J | A( 7 )|* 

We may specify an ensemble of Hermitian matrices by demanding that its joint probability 
density function is given by Qn. The nth correlation function of this ensemble is then 
defined to be R n : R™ — > [0, 00) where 

Rn(y) 1 / n N (yVj)dX N - n ( 7 ), (A.l) 

Z w {N-ny.J R N-n 

where y V 7 e R N is the vector formed by concatenating the vectors y € R™ and 7 € R ". 
By definition, R = 1. Here we take (A.l) as the definition of the nth correlation function; 
one can use the point process formalism to show that this definition is consistent with the 
definition derived in that manner. See [3] for details. 

We set fi n ,n' to be the n,n' entry of (Uq A )~ T , and we define q n := wq n . Using this 
notation we define the functions Sn,ISn and DSn '■ — ► R by 

2 Nb-l 2 Ar6_1 

SW(7.V) == ^ X] Qn{l)t\q~n'(i), IS N (-y,i') := - ^ /i n>n / e x q n {l) £\q~n>{i) 

n,n'— n,n'= 
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and 



Nb-l 



n,n'— 

The matrix kernel of our ensemble is then defined to be 



KnM) 



DSjvfr.V) ^(7,7') 

-s N tf,i) is N (~f,y) + 8(7, i). 



Theorem A.2. 



R n (y) = Pf [K N ( Yj , y y )] ; j, f = 1, 2, . . . , n. 
Our proof of this theorem begins by setting -q to be the measure on R given by 

T 



t=i 



where y\ , yi, • • • , Vt are real numbers and c\ , Ci , . . . , Ct are indeterminants and S is the prob- 
ability measure on R with point mass at 0. We will assume that T > N. As with Theorem 4 
in the main body of this paper, we will prove Theorem A.2 by expanding Z w ^ +r ^> jZ w in 
two different ways and then equating the coefficients of certain products of c\, c 2 , . . . , ct- 



Lemma A.3. 



= Pf (J + [VwK N (y u &,)]); t,t' = 1,2, ... ,T, 



where J is defined to be the 2T x 2T matrix consisting of 2 X 2 blocks given by 



J := 





' f 




St.e 


-1 





t,t=l,2,...,T. 



Lemma A.3 is proved in Section 8.1. 

For each JV>0we define 3„ to be the set of increasing functions from {1,2, ... ,n} into 
{1,2, ... ,T}. Given a vector y G R T and an element t 6 3^, we define the vector y t e R™ 

by y t = ij/t(i),yt(2), • • -,yt( n )}- 



Lemma A. 4. 



N ( T . 



(A.2) 



The proof of Lemma A. 4 is given in Section 8.2. 

Finally, we set K to be the 2T x 2T block matrix given by 

K := [Vc7c7K N (y t ,y t ,)]; t,t' — 1,2, ... ,T. 

From the formula for the Pfaffian of the sum of two antisymmetric matrices (see Lemma 7) 
and Lemma A.3 we have that 



Pf[J + K] = l + £ £PfK t , 

"=1 



(A.3) 
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where for each t <G 3^, K t is the 2n x 2n antisymmetric matrix given by 
K t = [y/Ci(j)Ct(j>)K N (y tU) , y t (j'))]; i'i' = 1,2,..., 5. 

Finally 

f n n 

j, j' = 1,2, . . . ,n, 



PfK t = ( n^-) jpf^^),^-,))]; 



and Theorem A. 2 follows from Lemma A. 4 by comparing coefficients of c\c 2 ■ ■ ■ c n in (A. 2) 
and (A.3). 



Appendix B The Pfaffian Cauchy-Binet Formula 

Theorem B.l (Rains). Suppose B and C are respectively 2 J x 2J and 2K x 2K antisym- 
metric matrices with non-zero Pfaffians. Then, given any 2 J x 2K matrix A, 

Pf (C- T - A T BA) Pf (B- T - ACA T ) 



Pf(C- T ) 



Pf(B" 



Proof. Let I 2 k and I 2 j be respectively the 2K x 2K and 2 J x 2 J identity matrices, and let 
O be the 2 J x 2K matrix whose entries are all 0. Then, an easy calculation shows that 



h j o " 




B T A" 




I 2 .7 


BA" 




A T B I 2K _ 




A T C- T _ 




o T 







B T 

o T 



O 

C T A T BA 



and similarly 



I 2 J 


-AC" 




B T A" 




hj O " 




B T ACA T O 


o T 


I2K 




A T C- T _ 




CA T I 2X _ 




o T C T 



Now, if D and E are 2N x 27V matrices and D is antisymmetric, then it is well known that 
Pf (EDE T ) = Pf D • detE. From which we conclude that 



Pf 



B T O 
O t C T A T BA 



Pf 



B ACA 

o T 



o 



The theorem follows since the Pfaffian of the direct sum of two even rank antisymmetric 
matrices is the product of the Pfaffians of the two matrices. That is 

Pf(B~ T ) Pf(CT T - A T BA) = Pf(B~ T - ACA T ) Pf(CT T ). □ 
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